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Although about a dozen high redshift {z larger than 2) starburst galaxies have 
been recently detected in the CO lines, spectroscopic detections of molecular gas 
of very young galaxies are still very difficult. The CO lines are usually optically 
thick, which limits greatly the increase of emission with redshift, as observed for 
the dust continuum. However the future instruments (LMT, ALMA, etc..) will 
allow large progress in this domain, and perspectives are discussed. Computations 
are based on a simple extrapolation of what is known of starbursting galaxies at 
lower redshift. 

1 Introduction 

One of the breakthrough due to recent progress on faint galaxies has been the 
inventory of the amount of star formation at every epoch (e.g. Madau et al. 
1996). The comoving star formation rate appears to increase hke (1 + z)* from 
2: = to z = 1, and then decreases again to the same present value down to 

2 = 5. But this relies on the optical studies, i.e. on the UV-determined star 
forming rates in the rest-frame. If early starbursts are dusty, this decrease 
could be changed into a plateau (Guiderdoni et al. 1997; Blain et al. 1999). To 
tackle this problem, independent information must be combined in a coherent 
picture, such as that coming from the far-infrared and millimeter domains: 
the cosmic IR and submm background radiation discovered by COBE (e.g. 
Puget et al. 1996; Hauser et al. 1998) yields an insight on the global past 
star-formation of the Universe, and the sources discovered at high redshift in 
the millimeter continuum and lines yield information on the structure of the 
past starbursts (Small et al. 1997; Hughes et al. 1998; Barger et al. 1999; 
Guilloteau et al. 1999). 

This review focus on the molecular content of galaxies, that we can de- 
duce from CO lines. Although the submillimeter continuum is more easy to 
detect at high redshift, it cannot without ambiguity trace the evolution of 
star formation as a function of redshift, because of identification problems, 
unknown redshifts, and unknown contribution of AGN. After summarising 
the present state of knowledge concerning CO emission lines, perspectives are 
drawn concerning the future surveys that will be conducted with the next 
generation of millimeter instruments. 
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2 CO Detections at High Redshift 



The detection of highly redshifted miUimeter CO hnes in the hyperluminous 
object IRAS 10214+4724 {z = 2.28, Brown & Vanden Bout 1992, Solomon et 
al. 1992), has opened this new field, to explore the history of star formation, 
and its efficiency. Although the first enormous derived H2 mass has now been 
revised (the source is considerably amplified by lensing), it is still surprising to 
find such huge amounts of CO molecules, especially since the gas is expected 
to have lower metallicity at high z. But theoretical calculations have shown 
that in a violent star burst, the metallicity could reach solar values very quickly 
(Elbaz et al. 1992). 

Today, more than a dozen objects have been detected in CO lines at high 
z: they are often gravitationally amplified, either being multi-imaged by a 
strong lens, like the Clovcrlcaf quasar H 1413+117 at z = 2.558 (Barvainis et 
al. 1994), the lensed radiogalaxy MG0414+0534 at z = 2.639 (Barvainis et 
al. 1998), or the possibly magnified object BR1202-0725 at z = 4.69 (Ohta 
ct al. 1996, Omont et al. 1996): or they are more weakly amplified by a fore- 
ground galaxy cluster, like the submillimcter-selcctcd hyperluminous galaxies 
SMM02399-0136 at z = 2.808 (Frayer ct al. 1998), and SMM 02399-0134at 
1.062 (Kneib et al. 2000). Often several high- J CO lines are detected, reveal- 
ing the high temperature of the gas (typical of a starburst at 60K). Higher 
temperatures are rare, as in the magnified BAL quasar APM08279+5255, at 
z = 3.911, where the gas temperature derived from the CO lines is ^ 200K, 
maybe excited by the quasar (Downes et al. 1999). Recently Scoville et al. 
(1997b) reported the detection of the first non-lensed object at z = 2.394, 
the weak radio galaxy 53W002, and Guilloteau et al. (1997) the radio-quiet 
quasar BRI 1335-0417, at z = 4.407, which has no direct indication of lens- 
ing. If the non-amplification is confirmed, these objects would contain the 
largest molecular contents known (8 — 10 • 10^° Mq with a standard CO/H2 
conversion ratio, and even more if the metallicity is low). The derived molec- 
ular masses are so high that II2 would constitute bc;tween 30 to 80% of the 
total dynamical mass (cf 4C 60.07, Papadopoulos et al. 2000), if the standard 
CO/H2 conversion ratio was adopted. The application of this conversion ratio 
is however doubtful, and it is possible that the involved II2 masses are 3-4 
times lower (Solomon et al. 1997). 

To search for primeval galaxies, already Elston et al. (1988) had identified 
extremely red objects that are conspicuous only in the near-infrared (R - K > 
5). Maybe 10% of the submm sources could be EROs (Small et al. 1999). A 
proto-typical ERO at z=lAA (Dey et al. 1999) has been detected in submm 
continuum (Cimatti et al. 1998), and in CO lines (Andreani et al. 2000). 
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Figure 1. H2 masses for the CO-detected objects at high redshift (full dots), compared 
to the ultra-luminous-IR sample of Solomon et al. (1997, open pentagons), to the Coma 
supercluster sample from Casoli et al (1996, filled triangles), and to the quasar 3c48, marked 
as a star (Scoville et al. 1993, Wink et al. 1997). The lines delineated by various symbols 
indicate the Icr detection limit at the IRAM-30m telescope of S(CO) = 1.0 Jy km/s with 
the 3mm receiver (hyphens), 2.0 Jy km/s with the 2mm (triangles) and 1.3mm (crosses) 
receivers, rms that we can reach in 6h integration. Note the absence of detected objects 
between 0.36 and 1 in redshift. The points at high 2 can be detected well below the 
sensitivity limit, since they arc gravitationally amplified. 



The CO line detections at high z up to now are summarized in Table |l|, 
and the molecular masses as a function of redshift are displayed in Fig. |l|. It is 
clear from this figure that only the very strongest objects are detected, and in 
particular the gravitationally amplified ones, but this will rapidly change with 
the new millimeter instruments planned over the world (the Green-Bank-lOOm 
of NRAO, the LMT-50m of UMass-INAOE, the ALMA (Europe/USA) and 
the LMSA (Japan) interferometers). It is therefore interesting to predict with 
simple models the detection capabilities, as a function of redshift, metallicity 
or physical conditions in the high-z objects. 

3 Perspectives with Future mm Instruments 

If so many sources have been detected now in the submillimeter domain at 
high redshifts, it is because around 1mm, starbursts have a larger apparent 
flux at z = 2 or 3 than z = 1 . The line emission does not have such a negative 
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K-corrcction, since in the low frequency domain, the flux of the successive 
hnes increases roughly as (optically thick domain), instead of i/'* for the 
continuum. Nevertheless the line emission is essential to study the nature 
of the object (the AGN-starburst connection for instance), and deduce more 
physics (kinematics, abundances, excitation, etc.). Given the gas and dust 
temperatures, the maximum flux is always reached at much lower frequencies 
than in the continuum, since the lines always reflect the energy difference 
between two levels; this is an advantage, given the largest atmospheric opacity 
at high frequencies. 



Table 1. CO data for high redshift objects 



Source 


z 


CO 
line 


S 
mJy 


AV 
km/s 


MH2 
10^° M© 


Ref 


SMM 02399-0134 


1.062 


2-1 


3 


500 


2* 


1 


0957+561 


1.414 


2-1 


3 


440 


0.4* 


2 


HRIO 


1.439 


2-1 


4 


400 


7 


3 


F10214+4724 


2.285 


3-2 


18 


230 


2* 


4 


53W002 


2.394 


3-2 


3 


540 


7 


5 


H 1413+117 


2.558 


3-2 


23 


330 


2-6 * 


6 


SMM 14011+0252 


2.565 


3-2 


13 


200 


5* 


7 


MG 0414+0534 


2.639 


3-2 


4 


580 


5* 


8 


SMM 02399-0136 


2.808 


3-2 


4 


710 


8* 


9 


6C1909+722 


3.532 


4-3 


2 


530 


4.5 


10 


4C60.07 


3.791 


4-3 


1.7 


1000 


8 


10 


APM 08279+5255 


3.911 


4-3 


6 


400 


0.3* 


11 


BR 1335-0414 


4.407 


5-4 


7 


420 


10 


12 


BR 0952-0115 


4.434 


5-4 


4 


230 


0.3* 


13 


BR 1202-0725 


4.690 


5-4 


8 


320 


10 


14 



* corrected for magnification, when estimated 

Masses have been rcscalcd to Ho = 75km/s/Mpc. When multiple images are 
resolved, the flux corresponds to their sum 

(1) Kneib et al. (2000); (2) Planesas et al. (1999) (3) Andreani et al. (2000); 
(4) Solomon et al. (1992), Downes et al. (1995); (5) Scovillc et al. (1997b); (6) 
Barvainis et al. (1994); (7) Prayer et al. (1999); (8) Barvainis et al. (1998); 
(9) Prayer et al. (1998); (10) Papadopoulos et al. (2000); (11) Downes et 
al. (1999); (12) GuiUoteau et al. (1997); (13) Guilloteau et al. (1999); (14) 
Omont et al. (1996) 
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3.1 Predicted Line and Continuum Fluxes 



To model high-redshift starburst objects, let us extrapolate the properties of 
more local ones: the active region is generally confined to a compact nuclear 
disk, sub-kpc in size (Scoville et al. f997a), Solomon et al. 1990, 97). The gas 
is much denser here than in average over a normal galaxy, of the order of 10'' 
cm~'^, with clumps at least of 10^ cm~'^ to explain the data on high density 
tracers (HCN, CS..). To schematize, the ISM maybe modelled by two density 
and temperature components, at 30 and 90K (cf. Combes et al. 1999). The 
total molecular mass considered will be 6 • 10^° Mq and the average column 
density N(H2) of 10^** cm~^, typical of the Orion cloud center. 

Going towards high redshift (z > 9) , the temperature of the cosmic back- 
ground Tbg becomes of the same order as the interstellar dust temperature, 
and the excitation of the gas by the background radiation competes with that 
of gas collisions. It might then appear easier to detect the lines (Silk & Spaans 
1997), but this is not the case when every effect is taken into account. To have 
an idea of the increase of the dust temperature with z, the simplest assumption 
is to consider the same heating power due to the starburst. At a stationary 
state, the dust must then radiate the same energy in the far-infrared that it 
receives from the stars, and this is proportional to the quantity T^^^^ — Tj^g, 
if the dust is optically thin, and its opacity varies in i^^, with /3 = 2. Keeping 
this quantity constant means that the energy re-radiated by the dust, pro- 
portional to T^ugt, is always equal to the energy it received from the cosmic 
background, proportional to T|^g, plus the constant energy flux coming from 
the stars. Since (3 can also be equal to 1 or 1.5, or the dust be optically thick, 
we have also considered the possibility of keeping T^^^^ — T^^^ constant; this 
does not change fundamentally the results. 

Computing the populations of the CO rotational levels with an LVG code, 
and in the case of the two component models described earlier, the predictions 
of the line and continuum intensities as a function of redshift and frequencies 
are plotted in Fig. ^. 

3.2 Source Counts 

To predict the number of sources that will become available with the future 
sensitivity, let us adopt a simple model of starburst formation, in the frame of 
the hierarchical theory of galaxy formation. The cosmology adopted here is 
an Einstein-dc Sitter model, = 1, with no cosmological constant, and Hq = 
75 km/s/Mpc, qo = 0.5. The number of mergers as a function of redshifts can 
be easily computed through the Press-Schechter formalism (Press & Schechter 
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Figure 2. Expected flux for the two-component cloud model, for various redshifts z = 0.1, 
1, 2, 3, 5, 10, 20, 30, and go = 0.5. Left arc the CO lines, materialised each by a circle 
(they are joined by a line only to guide the eye). Right is the continuum emission from 
dust. It has been assumed here that T'J'^^j. — T^^ is conserved (from Combes et al. 1999). 



1974), assuming self-similarity for the probability of dark halos merging, and 
an efficiency of mergers in terms of star-formation peaking at z 2 (i.e. 
Blain & Longair 1993). Also, the integration over all redshifts of the flux of all 
sources should agree with the cosmic infrared background detected by COBE. 
These contribute to reduce the number of free parameters of the modelling. 
To fit the source counts, however, another parameter must be introduced 
which measures the rate of energy released in a merger (or the life-time of 
the event): this rate must increase strongly with redshift (cf. Blain et al. 
1999). Once the counts arc made compatible with the submm observations, 
the model indicates what must be the contributions of the various redshift 
classes to the present counts. It is interesting to note that the intermediate 
redshifts dominate the continuum source counts (2 < z < 5), if we allow the 
star formation to begin before z = 6. At higher dust temperature, the counts 
are dominated by the highest redshifts (z > 5). 

Once these fits are obtained for the continuum sources, it is possible to 
derive also the counts for the CO line emission. The spectral energy distri- 
bution is now obtained with a comb-like function, representing the rotational 
ladder, convolved with a Planck distribution of temperature equal to the dust 
temperature, assuming the lines optically thick. The frequency filling factor 
is then proportional to the rotational number, and therefore to the redshift, 
for a given observed frequency. The width of the lines have been assumed to 
be 300km/s. The derived source numbers are shown in Fig. ^. Note how they 
are dominated by the high redshift sources. It is not useful to observe at A 
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Figure 3. Left Source counts for the CO lines at an observed frequency of 2rmn, assuming 
optically thick gas at Tex = 35K. The solid line is the total. Dash is the lowest redshifts 
(z < 2); Dot-dash, intermediate (2 < 2 < 5); Dots are the highest redshifts (z > 5). Right 
Same for A = 5mm. 



below Imm for high-^; protogalaxies, but instead to shift towards A = 1cm. 
4 Conclusion 

The search for CO hncs in high rcdshift galaxies is only beginning, there will 
be a real breakthrough in the next decade, when the sensitivity is increased 
by more than a factor 10 with the next generation millimeter instruments. 
CO lines are fundamental for our knowledge of star formation and efficiency 
at high z, to measure H2 masses, recognize AGN from starbursts, measure the 
dynamics and total masses of galaxies. 
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